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E.  I.  Jury^ 

SUMMARY 


In  this  study  a  simplified  analytic  test  of  stability  of  linear  discrete 

systems  is  obtained.  This  test  also  yields  the  necessary  and  sufficient 

conditions  for  a  real  polynomial  in  the  variable  z  to  have  all  its  roots 

inside  the  unit  circle  in  the  z-plane.  The  new  stability  constraints  require 

the  evaluation  of  only  half  the  number  of  Schur-Cohn  determinants^’  It 

is  shown  that  for  the  test  of  a  fourth-order  system  only  a  third  order 

determinant  is  required  and  for  the  fifth-order  only  two  determinants  are 

required.  The  test  is  applied  directly  in  the  z-plane  and  yields  the  minimum 

number  of  constraint  terms.  Stability  constraints  up  to  the  fifth-order  case 

are  obtained  and  for  the  n  order  case  are  formulated.  The  simplicity  of 

this  criterion  is  equivalent  to  that  of  the  Lienard-Chipart  criterion^  for  the 

continuous  case  which  has  a  decisive  advantage  over  the  Routh-Hurwitz 
4  5 

criterion  ’  . 


INTRODUCTION 

It  is  known  that  linear  time-invariant  discrete  systems  can  be 
described  by  constant  coefficient  linear  difference  equations.  These 
equations  can  be  easily  transformed  into  the  function  of  the  complex 
variable  z  by  the  z-transform  method.  One  of  the  problems  in  the  analysis 
of  such  systems  is  the  test  for  stability,  i.  e.  ,  to  determine  the  necessary 
and  sufficient  conditions  for  the  roots  of  the  system  characteristic  equa¬ 
tion  to  lie  inside  the  unit  circle  in  the  z-plane.  These  stability  tests 
involve  both  graphical  procedure  such  as  Nyquist  locus,  Bode  diagrams, 
and  the  root-locus,  and  analytical  methods  such  as  Schur-Cohn  or  Routh- 
Hurwitz  criteria.  Because  of  the  higher  order  determinants  to  be  evaluated 

II 

using  the  presented  form  of  the  Schur-Cohn  criterion,  many  authors  in 
the  past  have  used  either  a  unit  shifting  transformation^*  or  bilinear 
transformation  .  The  latter  transformation  maps  the  inside  of  the  unit 
circle  in  the  z  =  e  plane  into  the  left  half  of  the  w-plane  and  then 
applies  the  Routh-Hurwitz  criterion.  This  transformation  involves  two 


This  transformation  uses  p  =  z-1 
**  This  transformation  uses 


z  = 


w  +  1 
W  -  1 
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difficulties:  a)  algebraic  manipulation  for  higher-order  systems  becomes 
complicated,  and  b)  the  final  constraints  on  the  coefficients  in  the  z-plane 
become  unwieldy  and  require  algebraic  red>ctions  to  yield  the  minimum 
number  of  terms.  Because  of  these  limitations  this  criterion  is  not 
usually  used  for  systems  higher  than  second-order. 

A  recent  investigation  by  this  author  has  shown  that  the  evaluation 
of  the  Schur-Cohn  determinants  can  be  simplified  considerably  by  making 
use  of  the  real  coefficients  of  the  polynomial  in  z,  so  that  the  manipulation 
involved  in  testing  for  the  zeros  of  a  polynomial  are  comparable  to  those 
using  the  "transformed  (or  modified)  Routh-Hurwitz  criterion",  thus 
avoiding  the  bilinear  transformation 

The  study  in  this  paper  repre  sents  a  major  simplification  of  the 
earlier  work,  where  it  is  shown  that  only  half  the  number  of  determinants 
are  required  for  obtaining  the  stability  constraints.  This  simplification 
has  a  decisive  advantage  over  the  modified  Routh-Hurwitz  criterion  and, 
indeed,  higher-order  systems  can  easily  be  tackled  using  the  proposed 
stability  test. 


THEORETICAL  BACKGROUND 

In  this  section  we  review  the  simplifications  which  had  been  obtain¬ 
ed  in  an  earlier  publication  ^  and  explain  in  detail  the  manipulations 
involved. 

Schur-Cohn  criterion^’ 

If  for  the  polynomials 


F(z)  =  a^j 

+  a^z  +  a^z 
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are  different  from  zero,  then  F(z)  has  no  zeros  on  the  circle  |  z|  =1  and 
(J.  zeros  in  this  circle,  jji  being  the  number  of  variations  in  sign  in  the 
sequence  1,  Kl  .  I  I  ....  I  A^l  .  The  proof  of  the  above  theorem 
is  quite  involved  and  is  available  in  the  literature^’  This  criterion  was 
first  introduced  by  Cohn  in  1922,  and  since  that  time  neither  engineers 
nor  mathematicians  have  simplified  it  to  a  usable  form. 

For  a  system  of  order  n  to  be  stable,  all  the  n  zeros  of  its  charac¬ 
teristic  n  order  equation  must  lie  within  the  unit  circle,  i.  e.  ,  the 
sequences,  1,  |  Aj^  |  ,  |  A^  |  >  •  •  •  |a^|  must  have  n  variations  in  sign.  The 
stability  criterion  can,  therefore,  be  expressed  by  the  constraints^: 

1^1  <  0  ,  k  odd 

I  >  0  ,  k  even,  k  =  1,  2,  .  .  .  n  (3) 

For  a  discrete  or  a  sampled-data  system,  all  the  coefficients  of 
the  characteristic  equation  are  real.  Hence,  the  conjugate  sign  in  (2)  is 
superfluous.  It  is  the  utilzation  of  this  fact  that  leads  to  the  simplifica¬ 
tion  of  (2). 

As  noticed  from  (2),  the  highest-order  determinant  |  is  of 
order  2n,  while  the  characteristic  equation  is  of  order  n.  Hitherto  this 
constituted  one  of  the  discouraging  facts  in  widely  using  the  criterion  for 
higher-order  sampled-data  systems.  A  recourse  to  transformation  to 
other  planes  was  therefore  attempted  to  yield  easier  stability  tests. 

SIMPLIFICATION  OF  THE  STABILITY  CONSTRAINT  EQUATION®’ 

Since  all  a^^  in  equation  (2)  are  real,  the  matrix  can  be  written  as: 

where  the  superscript  T  denotes  transpose  and 

5  I'he  author  acknowledges  the  helpful  correspondence  with  Dr.  N.  H. 
Choksy  with  regard  to  the  material  in  this  section. 
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It  is  noticed  that  all  the  diagonal  terms  of  and  Qj^  are  equal  and  both 
arCf  symmetric  with  its  cross  diagonal.  It  is  this  characteristic  of  the 
Schur-Cohn  matrix  that  leads  to  the  following  simplification,  using  a 
unitary  transformation. 

Let  Ij^  be  the  k-order  identity  matrix,  I^^"^  the  k-order  permutation 

matrix, 


Uj^  the  2k-order  unitary  matrix 


'  (7) 


-1  T 

Let  Uj^  ,  then  by  actual  substitution  for  A^  ,  it  is 

readily  evident  that; 


(8) 
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One  can  easily  verify  that  P^=  I^Pj^  and  Qj^  =  ij^  q'^ 


(11) 


(12) 


(IS)"*^ 

(14) 

(15) 

(16) 


Thus  the  Schur=Cohn  determinant  |  Aj^j  is  reduced  to  the  product 
of  two  k-order  determinants  ’  ^  which  is  considerably  easier  to  evaluate 
than  the  direct  evaluation  of  the  2k-order  determinant  |  Aj^|  .  If  aj^  are 
complex,  then  this  simplification  is  no  longer  possible^’ 

THE  SYMMETRICAL.  PROPERTIES. OF  |  Xj^  I  ^  “  Y  |  ® 

Now  I  is  a  homogeneous  polynomial  of  dimension  k  in 

the  variables  a,  . . . ,  a^.  The  polynomial  |  -  Yj^|  is  identical  to  the 

polynomial  |  Xj^  +  Yj^  |  except  for  a  change  of  sign  of  those  monomial 
terms  which  have  an  odd  number  of  elements  from  Y  ,  i.  e. , 

^kf  =-^k  +  ®k 

■  '^kl  =^k  ■  ®k  (18) 

where  is  the  sum  of  all  monomial  terms  which  do  not  change 

(do  change)  sign  when  Yj^  is  replaced  by  -  Yj^  in  I  Xj^  +  Yj^  |  . 

identification  of  Aj^  and  (which  we  designate  as  the  stability 
constants)®,  ^ 

1)  Let  all  the  a.'s  in  the  matrix  Yj^  in  (12)  be  denoted  by  b.'s; 
then  expand  the  determinant  |  X,  +  Y,  I  in  terms  of  a.  and  b  . 

K  *  I  i 

2)  After  expansion,  examine  every  term  which  is  a  product  of 
a^'s  and  b.'s;  if  it  contains  an  even  number  (including  zero)  of  b.'s, 
then  it  is  assigned  to  Aj^;  otherwise  assign  the  term  to  Bj^. 

3)  After  collecting  the  terms  of  Aj^  and  Bj^,  replace  all  the 

b.'s  by  the  a. 's.  Hence 
1  1 


+  i  rom  (11)  and  (12)  by  repUcing  all  the  a's  of  Y  by  b's,  we  obtain  A,  and 
“Y  first  expanding  the  following  determinant.  1^ 
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and  for  stability  from  (3)  this  reduces  to 


I  I  ^  I  I  *  k  even 

^  •^k  I  I  ®k  I »  ^  k  =  .1,  2,  . . .  n  (20) 

Therefore,  the  application  of  the  Sch{ir.Cohn  criterion  now  reduces 
to  the  evaluation  of  determinants  up  to  order  n  only  for  the  n‘^  order 
polynomial. 

It  should  be  noted  that  the  last  simplification  is  useful  only  for 
design  procedures  where  the  coefficients  of  F(n)  are  given  in  other  than 
numerical  values.  However,  if  a^'s  are  given  in  numerical  values,  the 
use  of  equations  (17  and  18)  is  preferred. 

EQUIVALENCE  OF  THE  UAST.  CONSTRAINT  |  A J  4 1  B  |  TO  A  SIMPLE 

auxiliary  constraint^:  " 

The  constraint  which  we  will  introduce  involves  the  exclusion  of 
certain  real  roots  outside  of  the  unit  circle.  It  constitutes  a  necessary 
(but  not  sufficient)  condition  for  the  roots  of  F(z)  to  lie  inside  the  unit 
circle.  This  constraint  is  given  by 

F(z)  I  >  0  (1)  and  F(2)  I  >  0  for  n  even  (2a) 

^  ^  z=-l  <  0  for  n  odd  (2b) 

The  alternate  constrain  F(2)  |  <  0  and  F(z) .  <  0  for  n  even;  is 

2=1  0  for  n  odd 

also  possible.  However,  we  may  exclude  this,  without  loss  of  generality, 
by  always  letting  a^  be  positive,  in  which  case  to  satisfy  (1)  and  (2a,  2b) 
requires  that  F(2)  |  >  0  and  F(2)  |  >  0  for  n  even 

z=-oo  <  0  for  n  odd. 

Lemma  1.  U  (1)  and  (2b)  are  satisfied,  then  there  exists  at  least 
one  real  root  of  F(2)  =  0  between  plus  and  minus  one.  Also  the  total 
number  of  such  roots  is  odcf' 

Lemma  2.  If  (1)  and  (2a)  are  satisfied,  then  the  total  number  of 
real  roots  that  lies  between  plus  and  minus  one  is  zero  or  even^^. 
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To  show  the  equivalence  of  the  constraint  |  j  ^  |  B  )  to  the  above 
auxiliary  constraint  it  is  simple  to  distinguish  between  two  cases: 

(a)  n  is  odd:  Suppose  we  satisfy  the  constraint  constants  up  to 
\-l  ®n-l’  ^  generalization  by  Marden^of  the  Schllr-Cohn 

criterion  indicates  that  there  exist  (n-1)  roots  inside  the  unit  circle. 

The  arrangement  of  these  (n-1)  roots  (even  in  number)  inside  the  unit 
circle  is  one  of  two  alternatives.  (1)  The  first  alternative  is  that,  because 
complex  roots  appear  in  conjugate,  the  total  number  of  real  roots  between 
plus  and  minus  one  is  either  zero  or  even.  Now  if  we  impose  the  auxiliary 
constraint  (1)  and  (2b)  on  F(z)  we  find  that  the  last  single  real  root  from 
the  constraint  |  A^  )  <  (  j  should  lie  inside  the  unit  circle  from  Lemma  1. 
(2)  The  second  alternative  is  when  the  auxiliary  constraint  is  satisfied  in 
addition  to  the  first  (n-1)  constraints^  then  the  number  of  real  roots 
between  plus  and  minus  one  is  either  oi^  or  o^,  and  thus  in  this  arrange¬ 
ment  there  exists  a  single  complex  root  inside  the  unit  circle.  Since 
complex  roots  appear  in  conjugate,  the  last  constraint  |A  |<|b  |  is 
necessarily  satisfied.  Similarly  if  the  auxiliary  constraint  is  not  satisfied^ 
then  this  indicates  a  single  real  root  outside  the  unit  circle  and  thus  the 
last  constraint  is  also  not  satisfied. 

For  the  case  where  )  A^|  =  |  B^  |  this  indicates  a  real  root  on  the 
unit  circle  which  is  also  the  condition  of  the  auxiliary  constraint  when 
written  in  absolute  values  equated  to  zero.  Therefore,  we  have  shown  for 
n  o^  that  the  auxiliary  constraint  is  equivalent  to  the  last  constraint. 

(b)  n  is  even:  Suppose  we  satisfy  the  constraint  constants  up  to 
\-l  ®n-r  indicates  that  there  exists  (n-1)  roots  inside 

khe  unit  circle.  The  arrangement  of  these  (n-1)  roots  (odd  in  number) 
inside  the  unit  circle  is  one  of  two  alternatives.  (1)  The  first  alterna¬ 
tive  IS  that,  because  complex  roots  appear  in  conjugate,  the  total  number 
of  real  roots  between  plus  and  minus  one  is  either  one  or  odd.  Now 
if  we  impose  the  auxiliary  constraint  (1)  and  (2a)  on  F(z)  we  find  that  the 
last  single  real  root  from  the  constraint  I  A^  |  >  |  B^  |  should  lie  inside  the 
unit  circle  from  Lemma  2.  (2)  The  second  alternative  is  where  the 

auxiliary  constraint  is  satisfied  in  addition  to  the  first  (n-1)  constraints, 
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then  the  number  of  real  roots  between  plus  and  minus  one  is  either  zero 
or  even,  and  thus  in  this  arrangement  there  exists  a  single  complex 
root  inside  the  unit  circle.  Since  complex  roots  appear  in  conjugate 
therefore  the  last  constraint  ®nl  necessarily  satisfied.  Similar¬ 

ly  if  the  auxiliary  constraint  is  not  satisfied  then  this  indicates  that  a 
single  real  root  lies  outside  the  unit  circle  and  thus  the  last  constraint 
is  also  not  satisfied. 

For  the  case  when  I  A  I  =  I  B  I  this  indicates  a  real  root  on  the 

'  n  ‘  n 

unit  circle  which  is  also  the  condition  of  the  auxiliary  constraint  (written 

in  absolute  values)  when  equated  to  zero.  Therefore,  we  have  shown  for 

n  even  that  the  auxiliary  constraint  is  also  equivalent  to  the  last  constraint. 

Therefore,  for  stability  test  it  can  be  concluded  that  the  first 

(n-l)  constraints  of  the  A's  and  B's  should  be  satisfied  and  the  auxiliary 

constraint  is  then  equivalent  to  the  last  constraint  I  A  I  ^  I  B  I  .  This 

'  n '  '  n ' 

equivalence  has  been  checked  for  the  examples  discussed  in  this  note. 
Furthermore  in  the  next  sections  this  equivalence  will  be  demonstrated 
mathematically. 

THE  MODIFIED  STABILITY  CRITERION^  . 

Combining  the  previous  discussions  we  can  restate  the  stability 
criterion  in  a  modified  form  as  follows: 

A  necessary  and  sufficient  condition  for  the  polynomial  F{z)  =  a^ 

2  k  n  ^ 

+  A^z  +  a^z  +  .  . .  +  aj^z  +  . . .  +  a^z  ,  to  have  all  its  roots  inside  the 

unit  circle  is  represented  by  the  constraints 

I  1  <  I  \  I  for  k  odd 

and 

I  Aj^  j  >  I  Bj^  I  for  k  even,  k  =  1,  2,  .  .  . ,  n-l  (22) 

and  by  the  following  auxiliary  constraint. 


or 


F(z)|  >0 

z=l 


and  F(z)|  >0  n  is  even 

z=-l  <  0  n  is  odd 


for  a  >  0 
n 


F(l)  . 


F(-l) 


>  0 
<  0 


n  even 
n  odd 


for  any  a 
’  n 


(22a) 

(23) 
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MODIFIED  SCHtlR-COHN  CRITERION^ 

F rom  the  above  consideration  we  can  usefully  modify  the  Schur- 
Cohn  criterion  as  follows^; 

If  for  the  polynomial  with  real  coefficients 

F(z)  =  ao  +  +  ^2^^  +  . . .  f  a^z'^  .  a^  >  0  (24) 

satisfying  the  auxiliary  constraint,  all  the  stability  constants  Aj^  and 
(k  =  1,  . .  .  ,  n-1)  are  not  equal,  then  F{z)  has  no  zeros  on  the  circle  (  z(  =1 
and  {|i  +  1)  zeros  inside  the  unit  circle  for  n  even  and  p  ofid  as  well  as  for 
n  odd  and  p  even,  (p  is  the  number  of  variations  of  inequality  sign  in 
the  stability  constraints  [1,  (A^,  B^),  . .  .  ,  {A^_^,  .  Furthe  rmore, 

when  n  and  p  are  even  and  when  n  and  p  are  odd  the  number  of  zeros 
inside  the  unit  circle  is  p. 


REDUCTION  IN  THE  NUMBER  OF  DETERMINANTS  FOR  OBTAINING  THE 
STABILITY  CONSTANTS  Aj^'s;  a»d -^'s 

In  this  section  we  will  show  that  for  stability  test  only  about  half 
the  number  of  determinants  for  obtaining  the  A^^'s  and  Bj^'s  are  required. 
This  important  reduction  is  based  on  certain  properties  that  exist  between 
the  Aj^'s  and  Bj^'s.  We  will  indicate  these  properties  first  and  then  showi 
how  they  can  be  used  for  this  major  simplification 


1) 


A 


*k-i^fc+i 


k  =  2,  3,4  , 


n-l 


(25) 


The  above  equivalence  is  established  by  expansion  for  the  first 
few  values  of  "k"  and  can  be  generalized  similarly  for  any  other  value  of 
k  up  to  n-l.  When  k  =  n-l,  the  above  equivalence  is  reduced  to: 


> 

< 


A  ,A  <  B  -B 
n-c  n  n-2  n 


(26) 


Equation  (25)  is  the  key  identity  for  reducing  the  determinant  of  Aj^  and 
®k’  noticed  that  by  forcing  certain  restrictions  on  the  A's  and 

B's  before  and  after  a  certain  k  we  can  dispense  with  Aj^^  ^  ^^is 

will  be  best  illustrated  by  the  few  examples  to  be  discussed. 


Bn  =  (a^  +  a3  +  a5  + 


..){A 


n-l 


®n-l) 


■  ®n-l>  • 


(27) 

(28) 


10 


The  above  identity  can  be  easily  verified  for  the  first  few  values  of  "n" 
and  can  be  generalized  for  any  "n"^.  The  innportance  of  this  property 
lies  in  the  mathematical  proof  of  the  previous  section  as  follows: 


a2  ^ 

n  n- 


(aQ'^3,2+a^+. . . )  <  (aj^+a^+a-g. .  •  ) 


(29) 


or 


>  0 


\  -  \  <  0=±F(1).  F(-l)  J  J  (30) 

which  verifies  the  equivalence  |  A^|  ^  |  B^|  to  the  auxiliary  constraint. 

^‘^n-l'‘’®n-l^  "  ■^n-2^^0''’^2''’^4'*'*  ' '  ^  ■®n-2^^l'‘'^3'‘‘^5'''‘  *  *  ^ 


n  >  3 


From  (26)  we  can  write 


(A  ,-B  ,)(A  ,+B  ,)  ^  ^A  ..A  -B 

n-1  n-1''  n-1  n-1'  <  0  n-2  n  n-2  n  <  0 

using  (2  7)  and  (28)  in  the  right  side  of  equation  (32) 

A  .,(a  +a  +a.+a.. ..  )(A  ,-B 


(31) 


(32) 


>0 


*n-2'“0'“2'“4“^6‘*‘ ''“n-1  "n-l 


^“®n-2K'^^3'^^5'  •  •  ^^^n-r®n-l^<0 


or 


(33) 


Comparing  equation  (34)  with  the  left  side  of  equation  (32)  we  obtain  the 
following  identitv: 


>  0 


^^n-l'^®n-l*  <  0— ^n-2^^0‘^^2'^^4'*''  '  '  ^■®n-2^ V^3'*'^5*  *  ‘  ^  <  0  ’  "  -  ^ 

(35) 

The  above  relationship  can  also  be  derived  directly  for  the  first 
few  values  of  "n"  and  can  be  generalized  for  any  "n".  Furthermore 
the  above  property  can  also  be  shown  as  an  equivalence  relationship 
which  can  be  written  as  a  recurrence  equation. 

^  Vl'^®n-1>  =  ^n-2^^0^^2'^^4-  *  *  )-®n-2^ V^3'^^5*  ’  ’  ^ 

The  use  of  this  property  lies  in  the  fact  that  one  can  obtain  directly 

(A  ,+B  )  from  the  previously  obtained  A  _  and  B  and  also  to 

verify  the  first  property  for  the  upper  limit,  i.  e. ,  k=n~l  (See  Appendix). 

The  above  three  properties  in  combination  with  the  preceding  discussion 
will  now  be  used  in  obtaining  the  new  stability  constraints  for  low  order 
systems  and  then,  by  generalization,  to  obtain  the  constraint  for  an  n- 
order  system.  * 

+  I'o  establish  the  identity  for  the  general  case,  it  is  easier  to  show 
the  following  equivalent  conditions  for  (2  7)  and  (28): 

A  +B  =  F{1)(A  ,-B  ,)5  A  -B  =  F(-1){A  ,-B  ,) 

n  n  n-1  n-l'  n  n  \  ir  ^-1  n-1' 
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Examples  of  Low  Order  Systems: 

We  will  apply  the  reduction  properties  wherever  they  are  appli¬ 
cable  to  n  =  2,  3, 4i  5,  and  then  obtain  the  stability  constraint  for  any  n. 

We  will  assume  a^  >  0.  This  can  be  easily  done  by  multiplying  F(z) 
by  minus  one,  if  necessary. 

(a)  n  =  2,  F(z)  =  ajj  +  ajZ  +  a^z^,  a^  >  0  (37) 

The  stability  constraints,  using  equations  (22  and  22a)  are  given 

^^2  ®ll  ■  ■  (38) 

VV^2^°*  VV^2^°'  >182!  (39) 

One  could  also  remove  the  absolute  sign  from  equation  (38)  for 
if  a^  is  negative  with  magnitude  larger  than  a^,  equation  (39)  will  not  be 
satisfied.  However,  we  may  leave  the  absolute  sign  in  order  to  discontinue 
the  stability  test  if  (1)  is  violated.  , 

(b)  n  =  2  ,  F(z)  =  +  aj*  +  a^z^  +  a^z^  ,  a^  >  0  (4O) 

The  stability  constraints  from  the  modified  stability  equations  (22)  and 
(22a)  are  given  by  the  following  inequalities: 


1^0  1  ^  ^3  *  1 

l^il  <  iBjl 

(^) 

-  ^3^1 

^  1  ^0*2  ”  ^1*3  ^  ’  1  ^Z  ^  ^  1  ®2  ^ 

(42)'^ 

F(z)  1  >  0  , 

Z  =  1 

F(z)(  <  0  ,  or  1  A  1  <1  B  1 

z=l  ^  ^ 

(43) 

To  show  how  to  obtain  and  for  this  case: 


1)  Expand  the  determinant  of  Eq.  (20a)  for  k=2  and  n=3,  as  follows: 
|X^  +  Y  I  = 


“0  +  '>2 


a,4  bj 


'  “0^  +  “o'’2  -  »l'‘3  -  '>3^ 


3  “0 

2)  To  identify  and  B^,  follow  procedure  (2)  on  page  6  ,  to  obtain 

A  -  2,2 

■  ^3 

®2  ^  ®0*^2  "  *1*^3 

3)  Replace  all  the  b’s  of  A  and  B  by  the  a's  to  get: 

.22^ 

^2  *0  "‘“3 

®2  =  V2  -  »1*3 

A  similar  procedure  is  used  for  obtaining  the  A,  's  and  B,  '(s)  for 
any  "k"  and  "n".  ^  ^ 
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Reduction  of  the  constraint  conditions: 

From  the  first  property  of  the  Aj^  and  ,  we  may  write  in  t^’  'ase 

(44) 


Since  Bj  is  positive  (i.e.,  a^  >  0) ,  we  may  write  for  the  right  side 

■^1 

-g—  A^  >  ®3  • 

Using  conditions  (41)  and  (43)  in  combination  with  {42a)  the  new  stability 
constraints  are: 

B^  <  0 

U3I  <  I  B^l  ,  or  F(l)  >  0,  F{-1)  <  0 
From  the  second  property  in  equation  (28) 

B3  =  (aj  +  -  Bp 

Since  +  a^  is  positive  (from  Eq.  4  7)  because  it  is  equal  to 

F(l)  -  F(-l)  >  0  ,  therefore  B^  is  negative  only  when  -^2  " 

the  third  property  of  the  stability  constants  it  is  readily  seen  that 

A2  +  B2  <  0  is  identically  satisfied  from  the  first  and  third  conditions. 

Thus  the  simplified  form  of  the  stability  constraint  for  n*:3  reduces  to 

I%l<a3 


(42a) 

(45) 

(46) 

(47) 

(48) 


0 


^3  <^0^2-V3 


^0  ^  +  ^2  ^3  ^  ^0  ■  ^  ^2  ■  ^3  <  ° 

Stability  diagrams  for  a  second  and  third-order  case  are  presented  in 
Figures  1  and  2.  ^ 


(^)  H  ~  F{z)  -  a.  +  a  z  +  a_z^  +  a_z^  +  a  z"^ 
u  1  z  3  4 


a .  >  0 
4 


The  stability  constraint 

l^ol  l^ll  <  l^il 


(49) 

(50) 

(51) 

(52) 

(53) 


‘0 


-  a 


■|  >  ho 


^3  -  ^1^4 


A^l>  I 


^^^0^2^4  +  ^1^3^  -  ^0^4  -  ^2^4  " 


(54) 

3 


-^0^1^3 


(55) 
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or 


(55) 


4 


■3^ 


IA3I  <  IB3I 


F(z),  >  0  ,  F(z)  >  0  ,  or  I  A  I  >  I  B. 

'z=l  U-l  4'  '4 


(56) 


Reduction  of  the  constraint  equations 

Using  property  (1),  we  may  write  as  in  the  previous  case: 

•^2  ^  ^  "5^  -^3  ^  -  (57) 

Equations  (53),  (54),  and  (55)  are  now  equivalent  to  (53),  <  0  and 

(55).  Using  B3  <  0  with  (55),  we  finally  obtain  the  reduced  constraints: 


(58) 


A3  -  B3  >  0,  A3  +  B3  <  0  (59)+ 

F(l)  >  0,  F(.l)  >  0  (60) 

It  is  noticed  that  for  the  fourth-order^^  case  only  one  determinant 
of  third-order  for  obtaining  A«  and  B_  is  required.  All  other  conditions 
are  very  simple. 

(d)  n  =  ^,  F(z)  =  a^+a^z+a^z  +a3Z  +a^z’+a3Z^,  a^  >  0  r  (61) 


The  stability  constraints  in  symbolic  form 


1  Ajf  <  B^  ,  =  ag  >0 

(6  2) 

lA^I  >  IBJ 

(6  3) 

IA3I  <  IB3I 

(6  4) 

IA4!  >  |B^| 

(6  5) 

F(l)  >  0  .  F(-l)  <0,  or  1  A3I  <  1  B3I 

(66) 

Constraints  in  (59)  can  also  be  written  as  B3<0,  |  A  |  <  Ib- |  .  The  latter 
has  an  advantage  in  numerical  testing  if  B3<  0  is  vicflated.  Then  the  test 
could  be  discontinued  without  having  to  calculate  A3. 

++An  alternate,  form  which  is  more  advantageous  for  design  can  be  obtained 
for  the  fourth  order  case.  This  form  can  be  easily  obtained  by  using 
properties  (1)  and  (2).  It  is  given  as  follows:  (1)  A-  <0,  (2)  A,  <  -  I  B,|  , 

(3)  A3  -  B3  >  0,  (4)  F(l)  >  0,  F(-l)  >  0.  •  ^  2  '  2' 

In  this  case  only  one  third  order  equation  in  (3)  and  one  second  order 
equation  in  (2)  are  to  be  solved,  while  in  the  former  case  two  third  order 
equations  are  to  be  solved.  It  should  be  noted  that  when  (2)  is  satisfied, 
relationship  (1)  becomes  redundant. 
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Reduction  of  the  constraint  relationship;  ® 

We  may  keep  in  this  case  condition  (63)  but  we  eliminate  (64) 
by  using  the  first  property,  i.  e.  , 

a\  <  b] 


:  (64a) 


<  0,  <  0 


A4  -  >  0,  +  B^  >  0 


With  negative  because  of  (62),  the  stability  constraint  for  (62),  (64), 
(64a)  and  (65)  becomes:  (62),  (63),  A^  >  0,  and  A^  >  .  Using  A^  >  0  “ 

in  ( 65)  i.  e. ,  (A4-B^)(A^+B^)  >  0  we  obtain  A^  -  >  0  and  A^  +  B^  >  0. 

Furthermore,  since  A2  =  A^  -  B^  is  to  be  negative  from  (62),  it  is 
readily  satisfied  if  the  second  constraint  is  equivalent  to  A  -  B  <0, 

■^2  ®2  Finally,  we  obtain  for  the  reduced  form  the  following: 

(67) 

(68) 

F(l)  >  0,  F(-l)<  0  (69) 

It  is  noticed  that  in  this  case  only  the  second  and  fourth  order  determinants 
are  required.  One  may  also  use  a  different  form  of  reduction  by  eliminat¬ 
ing  (68).  However,  this  will  not  yield  much  simplification  over  the 
previous  form  because  a  fourth-order  determinant  with'a  third-order 
determinant  is  then  required. 

The  above  discussion  can  be  generalized  for  an  "n"  which 
finally  reduces  to  the  following  simplified  criterion. 

The  New  Stability  Criterion: 

A  necessary  and  sufficient  condition  for  the  polynomial  ' 

F(z)  =  Uq  +  a^z  +  a2Z^  +  . . .  +aj^z^+.  .  .  +a^z”  with  a^  >  0,  to  have  all  its 
roots  inside  the  unit  circle  is  represented  by  the  following  constraints 
for  n  even  and  n  odd  respectively: 


Note  A  =  a/-  B  ^ 
^11 
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+ 

n  even 

n  odd^ 

A^  -  B^  <0,  A^  +  <  0 

A^  “  >  0,  A^  +  B^  <  0 

A4  -  B^  >  0,  A^  +  >  0 

A^  -  Bg  <  0,  Ag  +  Bg  >  0 

^6- 

A^  -  B^  >  0,  A^  +  B^  <  0 

•  • 

•  « 

•  • 

^  _g  >0  for  n=4k 

n-1  n-1  <  0  for  any  other  n 

•  • 

A  -  B  ^ 

n-1  n-1  <  0  for  any  other  (n-1) 

^  ^  g  <0  for  n  =  4k 

n-1  n-1  >  0  for  any  other  n 

A  +B  ^  ®  n-l=4k 

n-1  n-1  <  0  for  any  other  (n-l) 

k  =  1,2,  3. .  . 

k  =  1,  2,  3. .  , 

F(l)  >  0,  F(-l)  >0 

F(l)  >  0,  F(-l)  <0 

Alternate  Forms: 

An  alternate  equivalent  method  which  is  of  advantage  if  the  stability 

constants  evaluation  is  carried  out  by  methods  other  than  a  computer  is 

hereby  presented. 

n  even 

n  odd 

|Aj^  I  <  Bj^  ,  Bj  =  a^  >  0 

1  -^2  1  ^  1  ^2  ^  '  ^2  ^  * 

IA3I  <  1.33 1  ,  B3<0 

1  A4I  >  1  A^  >  0 

1  A5I  <1  B3I  .  B3  >  0 

|A^|  <  1  B^l  ,  B^<  0 

•  • 

**  *  •• 

I  A  ,1  <|b  ,  I  ,  B  <0.  for  n=4k  .  I^Ir  |  a  >0,  for  n-l=4k 

''-I  n-l*  n-l>0,  for  any  other  '  n-1  <  I  ^n-l' ’  ^n-KO,  for  any  other 

"  n 


k=l,  2,  3.  ,  , 

F{1)  >  0,  F(-l)  >  0 


k=l,  2,3... 

F{1)  >  0,  F(-l)  <  0 


■^Note  that  A  +  B,  and  A, 
(17)  and  (18)^  ^ 


-  can  be  obtained  directly  from  equations 
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In  the  above  case  the  identification  of  |  Aj^|  and  |  Bj^|  from  the 
determinant  j  -  Yj^|  could  be  used.  This  criterion  can  be  applied 
for  design  purposes  when  the  of  F(z)  are  given  other  than  numerically. 
Illustrative  Examples 

To  illustrate  the  stability  test,  we  choose  two  problems,  one 
involving  design  and  the  other  a  numerical  test  of  a  polynomial. 

1.  The  design  problem  is  concerned  with  obtaining  the  maximum 

3 

allowable  value  of  k  (the  gain)  of  a  feedback  sampled-data  system^^ 
shown  in  Figure  3,  to  be  stable. 


The  overall  transfer  function  is  :given  as: 

C(z)  ^  G(z) 

“rR  -  r+  G|z) 

where 


(70) 


G(z)  =  z-transform  of  the  forward  path  transfer  function  G(s) 


e 


-0.  5s 


0.  5k 
s  +  0.  5 


0.118k 
-  0  .  882) 


(71) 


For  stability  we  have  to  examine  the  denominator  of  Eq.  (70)  i.  e. , 
1  +  G(z)  =  z^  -  0.  882  z^  +  0.118k 


1  +  G(z) 
(72) 


The  above  equation  is  a  third-order  polynomial  in  z,  and  thus  we  can 
apply  the  stability  tests  obtained  earlier  for  n  =  3.  In  this  case 

^3  “  ^2  ~  ~  ^0  ~  ^  positive. 

1)  la^l  <  a^,  0.118k<l,  k  <  8.  47  (73) 

2  2 

^0  ”  ^3  "  ^0^2  ^  yields  k  <  5.  75  or  k  ^  -13.  21  (for  positive 

feedback) 

^0  ^  ^1  ^  1  -0.  882  +  0.118  k  >  0  is’^tisfiedfor  any 

positive  k  >  0 

^0  ~  ^1  ^2  ■  ^3  +  0.118  k  <  0,  k  <  15.  9  (74) 
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Therefore,  the  maximum  allowable  gain  for  stability  is  the  lowest  value 

which  is  in  this  case  k  =5.  75. 

max 

2.  Test  for  stability  the  following  polynomial: 

F{z)  =  z^Zz^-O.  5z-0.  95  (75j 

The  above  polynomial  is  again  of  third  order,  i.  e.  ,  n  =  3,  thus  we  apply 
the  stability  constraints  for  this  case,  a-=l,  a  =  2,  a  =-0.  5,  a  =  -0.95 

1)  la^l  <  a^,  0.95<1 

2)  aj  -  a^  <  a^a^  -a^a^  ,  0.  95^  -  1  >  -1.  90  +  0.  5  (76) 

or,  -0. 1  >  -1.  40 

The  second  condition  is  violated,  thus  there  exists  at  least  one  root  out¬ 
side  the  unit  circle  and  thus  the  system  is  unstable.  To  determine  the 
number  of  roots  outside  the  unit  circle  from  the  modified  Schtlr-Cohn 
criterion,  (see  p.  10)  we  also  have  to  examine  the  sign  of  the  last  condition 

3)  a^+a^+a^+a^  >0,  -0.  95  -0.  5  +  2  +  1  >  0 

^0"V^2"^3  <  0*  “0-  95  +  0.  5  +  2  -  1  >  0  (77) 

The  last  condition  |  <  |  Bj|  is  violated.  Now  the  number  of 

changes  of  sign  of  1,  Aj^'s  and  Bj^'s  are  the  number  of  roots  inside 

the  unit  circle.  In  this  case,  the  sign  changes  are,  1,  A?  -  <  0. 

22  22  11* 

^2  ~  ^2^  "^3  "  ®3  ^  There  are  two  changes  of  sign,  and  since 

only  three  roots  exist,  therefore  only  a  single  real  root  exists  outside 

the  unit  circle. 


CONCLUSION 

From  the  preceding  discussion,  it  is  shown  first  that  in  the 

I! 

original  Schur-Cohn  or  the  modified  Routh-Hurwitz  criterion,  the  number 
of  determinants  required  for  the  stability  is  almost  halved.  The  use  of 


rt\ 


18 


the  criteria  for  both  design  of  discrete  systems  as  well  as  for  testing 
roots  of  a  polynomial  inside  the  unit  circle  is  illustrated.  This 
criterion  will  now  be  useful  in  many  applications  such  as  the  stability 
test  of  difference  equations  with  constant  and  periodically  varying 
coefficients,  in  nonlinear  discrete  systems  for  the  stability  study  of 
limit  cycles,  in  the  design  of  digital  computers,  in  the  stability  test 
of  linear  systems  with  randomly  varying  parameters  and  in  many 
other  applications.  Thus  it  is  hoped  that  this  criterion  will  find  many 
applications  in  various  fields  in  addition  to  the  above  and  its  use  by 
engineers,  physicists  and  mathematicians  will  be  greatly  enhanced. 
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FIG.  I  STABILITY  DIAGRAM  FOR  A  SECOND  ORDER  CASE 
F(z)=aoZ  +  Q|Z  *02^^  ,  Q2=I 


20a 


-  -CSV; 


ZERO-ORDER  PURE 

HOLD  DELAY  PROCESS 


T  =  0.25  sec 


« 


FIG,  3  A  SAMPLED  DATA  FEEDBACK  SYSTEM 


APPENDIX 


General  Proof  of  the  Properties  of  the  Stability  Constants  Aj^'s  and  Bj^'s 

In  this  appendix  the  second  and  the  third  properties  of  the  stability 
constantsti.  e.  >£qns*(2 7),  (28)  and  (31),  will  be  mathematically  proven 
and  from  these  two  properties,  the  limiting  case  of  the  first  property, i.  e. , 
Eqn,(26)  will  be  demonstrated.  A  heuristic  argument  will  be  presented 
to  indicate  that  relationship  (25)  is  valid  for  all  k  =  2,  3,  .  . . ,  n~l. 

Proof  of  the  Second  Property: 


The  second  property  as  indicated  in  Eqns.  (2  7)  and  (28)  is  given 
follows: 


as 


A„  =  (a„t  aj  t  +  aj,  t  .  . .  )(A„_j  -  B„_j) 

B„  =  (ai  +  a^t  ajt  a  ,  a  J.  2 

To  show  this  property  it  is  simpler  to  manipulate  the  following 
equivalent  relationship,  which  is  obtained  by  adding  and  subtracting 
the  above  two  equations. 


A  +  B  =  (a  +a,+a-. 
n  n  '  o  1  2 


=  <^-»l"-2-  •  •  •  •  -  B„.i) 


The  above  can  be  also  written  as: 


An  "  ^n  =  (A^_^-B^_^) 


(1) 


A  -  B  =  F(-1)(A  -B  ) 
n  n  '  ''  n-1  n-F 


(2) 


au  t-  acknowledges  the  aid  of  Mr.  Jean  Blanchard 

the  discussions  of  this  appendix. 
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We  will  first  demonstrate  relationship  (1)  and  following  the  same 
procedure  relationship  (2)  can  be  similarly  demonstrated.  The  proof 
consists  ofdeterminant  manipulation  and  in  particular  using  the  follow¬ 
ing  property. 

"The  value  of  the  determinant  is  unchanged  if  the  elements  of  any 
row  (column)  are  replaced  by  the  sums  or  (differences)  of  the  elements 
of  that  row  and  the  corresponding  elements  of  another  row  (column).  " 

To  show  relationship  (1),  we  write  first  the  determinant  A  +B 

/II  n  n 

as  follows: 


A  +B  =  X  +  Y  1  = 
n  n  '  n  n' 


a  +a,  a,+a . a  ,+a 

o  1  12  q-1  q 


1+1 


i+2 


a_  . a  _+a  , 

2  13  q-2  q-1 


a  ,+a 
n-1  n 


n-2 


V<+l%+i-2-***  Vf+1 


V<+1'^V<-1  ^n-f 


. Vf-1 


a  0  0 

n 


0  Oa 


(3) 


We  will  show  the  equivalence  of  (1)  for  the  general  row  and 
column  in  the  above  matrix,  by  concentrating  only  on  rows  f-1, 
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I,  1+1  and  columns  q-1,  q  and  q+1  as  follows: 


row  f  -1 

^q-l'‘‘\+f-3 

-f +1^  ^,q+Jt-2 

\-f+2'’’\+f-l 

row  f 

*q-f  -2 

a  -.,+a  .. 
q-f  +1  q+f 

row  1+1 

\-t  -1 

^q-1 -l^^q+f 

a  +a 
q-f  q+f+1 

column  q-1 

column  q 

column  q+1 

(4) 


with  a|i  =  0 
for  >  nor  )j.<  0 


Similarly  we  obtain  the  same  rows  and  columns  for  the  determinant 


n-1  n-1 


General  Coefficients  inMatrix  A  ,-B  , 
_ _ _  n-1  n-1 


^q-f+2 

row  f  -1 

— 

— 

— 

®q+f -2 

^q+l -1 

*q+f  -1 

Vf-1 

Vf+1 

row  f 

— 

— 

-1 

\+l 

^q+f+1 

Vf  -2 

Vf-1 

row  f+1 

Vf+1 

^q+f+2 

column  q-1 

column  q 

column  q+1 

(5) 
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To  demonstrate  (1)  we  perform  the  following  matrix  manipulations: 
(1)  In  the  determinant  A^+B^  ,  add  up  to  the  last  row,  all  the  previous 
rows.  We  then  notice  that  all  the  coefficients  of  the  last  row  after 
this  adding  become, 

p=n 

1 

p=o 


The  new  determinant  after  factorization  of  F(l)  becomes. 


n  columns 

The  elements  in  determinant  D  are  the  same  as  the  initial  determinant 

A  +B  . 
n  n 

(2)  In  (6)  we  subtract  the  colunms  2  from  1,  ...  the  column  q  from 

q-1,  .  .  .  and  the  column  n  from  n-1.  By  so  doing,  it  is  noticed  that 

except  for  the  coefficients  of  the  column  n  and  row  n  which  is 

equal  to  1,  all  other  coefficients  of  the  row  n  are  equal  to  zero. 

Therefore  the  matrix  in  A  +B  can  be  written  as: 

n  n 
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To  identify  D,  with  the  matrix  A  ,-B  ,,  in  (5),  we  have  to 

n-i  n-i 

eliminate  the  encircled  coefficients  in  the  determinant  D^. 

(3)  To  show  the  above  we  rewrite  for  simplcity  only  the  coefficients 
in  the  column  q  of  the  Dj  determinant  as  follows: 
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To  cancel  the  encircled  coefficients  in  or  in  above, we 
perform  the  following  operations:  Add  up  row  1  to  2,  we  cancel  out 
the  encircled  term  in  (2),  then  considering  the  new  row  2  obtained 
and  adding  this  row  to  row  3J*,  we  cancel  out  the  encircled  terms  in 
row  (3),  we  continue  this  process  to  cancel  all  the  encircled  terms. 
Finally  we  obtain  the  column  q  in  as  follows; 

Column  q  in  Dj^ 


row  1 

q-1  q+1 

row  2 

^q-2"V2 

row  3 

^q-3"V3 

row  < 

^q-l ’^q+l 

• 

• 

i  • 

• 

row  n-1 

0  ,  if  q  1,  when  q=l,  it  is  equal  to 

or  q  ^  n-1  it  is  equal  to 

a 

o 

Comparing  this  column  with  column  q  in  Eqn.  (5),  we  readily 

establish  the  equivalence  which  is  valid  for  all  q  columns  and  i 

rows:  Thus  the  identity  between  A  -B  ,  and  D,  is  established. 

n-i  n-1  1 

Therefore  Eqn.  (1)  is  verified. 

Following  the  same  procedure  with  the  appropriate  operations, 
Eqn.  (2)  can  be  similarly  verified  and  thus  the  second  property  is 
•  demonstrated. 
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Proof  of  the  Third  Property  of  and 


The  third  property  as  given  in  Eqn  (31)  is  represented  as  follows: 


A  ,+B  ,  =  (a  +a-+a  .+ 

n-1  n-1  o  2  4 


)  '^n-2"^V^3'^  •  •  •  ■'■^Zp+l'*'’  •  •  ^®n-2 


The  above  relationship  is  also  equivalent  to  the  following, 


In  this  discussion  a  rigorous  proof  of  the  above  relationship  will 
be  obtained  from  which  the  above  property  is  established.  Furthermore 
by  combining  the  third  property  with  the  second,  we  will  establish  the 
limiting  case  of  the  first  property,  i.  e.^  k=  n-1. 

The  proof  will  be  based  on  determinant  manipulations,  by  using  the 
same  properties  as  in  the  previous  case. 

The  determimant  of  A  ,,+B  ,  can  be  written  as: 

n+1  n-1 
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column  2q+l  column  2q+2 


a  +a-, 
o  2 

ai+as 

^2q'‘’^2q+2  ^2q+1^2q+3 

a  _+a 
n-2  n 

A  +B  = 

n-1  n-1 

“3 

^0+^4 

^2q-l''‘^2q+3 

• 

a  -1 
n-3 

|X  ,+  Y  ,1  = 

'  n-1  n-1' 

> 

^5 

^2q-2'^''2q+4 

a  A 

n-4 

row  2| 

^21+1 

^2q-2|+i*’^q+2l+2 

• 

• 

row  2|+1 

^21+2 

^2q-2|'''  ^2q+2l+2 

• 

• 

row  2|+2 

^21  +  3 

« 

^2q-2|-l''’^2q+2l  +  3 

• 

• 

a 

n 

0  0 

a 

o 

For  the  determinant  A  _  +  B  we  can  write  in  a  similar 

ri* b  • 

fashion,  however  we  concentrate  on  the  general  rows  and  columns 
as  shown: 


Aiatrix  of  A  _  +  B  ^ 
n-2  -  n-2 


with  aji  =  0  if  <  u  or  >  n 
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Procedure  for  the  Proof: 


1.  Using  the  matrix  A  +  B  ,  w 

3.  5.  7. .  .  I^2p+lJ‘ 

of  the  new  matrix  obtained  becomes: 


we  add  to  row  1,  the  rows 
readily  seen  that  the  first  row 


column  1 

column  2 

column  2q+l 

column  2q+2 

row  1 

y  2p+l 
- - 

Z“2p 

Z  ^2p+l 

with 


1  *2p  ”  ■  • 

2,*2ptl  = 


2.  Then  subtract  from  the  columns  2q+l  °  P  n=2p+l 

’  o^  q<  p-1  if  n?:  2p 

the  columns  2q-l.  Similarly  we  subtract  from  the  columns  2q, 
the  columns  2q-2,  this  operation  being  performed  step  by  step. 

For  instance  if  n=6,  we  first  subtract  the  column  3  from  the 
column  5.  Then  the  column  1  from  3.  Similarly  we  subtract  the 
column  4  from  6  and  then  the  coli^n  2  from  4.  3y  performing  this 
operation,  we  notice  that  the  first  row  contains  all  zeros  except  for 
the  first  and  second  columns,  where  the  coefficients  are  now 

^^2p  ^^2p+l  ’  noting  that 

I 

I 

5  [F(l)  +  F(-l)  ]  =  ^  ,  and  ^  [F(l)-F(-1) )  =  , 

the  first  row  of  A  +B  ,  becomes; 

n-i  n-i 


column  1 


column  2 


A  ,+B  * 

n-1  n-1 

Row  1 

i  [F(1)+F(-1)] 

i  [F(l)-F(-1)] 

000.  .  .  000.  .  .  00 

4 

3.  Add  column  2  to  column  1,  then  in  the  new  determinant  multiply 

column  2  by  Z,  and  divide  the  determinant  by  2,  and  then  subtract 

column  1  from  colximn  2,  we  get  for  the  determinant  A  .+B 

n+1  n-1 


column  1 


column  2 


in  the  determinant  A,  the  following  coefficients  in  the  general  row 
and  column  appear: 

column  2q+l  column  2q+2 


Roy/  2< 

^2q-2<+l'^^2q+2<+l 

"^2q-2l-r^2q+2f-l 

^2q-2f  +  2'‘’^2q+2f  +  2 

■^2q-2f ■^2q+2f 

Kl  < 

o<q< 

^2q-2l'*'^2q+2f +2 

^2q-2<+l'^^2q+2f  +  3 

Row  21+1 

^2q-2<-2  ^2q+2f 

■^2q-2f -r^2q+2l  +  l 

We  expand  the  previous  determinant  with  respect  to  the  first 
row  to  obtain; 
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where  and  Dj  are  the  appropriate  determinant  obtained  from  i . 

If  we  can  show  that  and  D,  =  A.  ,+B.  ,  ,  then 

we  complete  the  proof  of  the  third  property. 

4. 

and  Dj,  as  follows:  is  obtained  from  A  as. 


We  will  demonstrate  first  the  equivalence  between  A  +B 

n-2  n-2 


column  1 

column  2 

Row  1 

(a^+a3)+a4 

a,+a_-a- 
1  0  3 

Row  2 

a  +a>-a^ 
o  6  4 

column  2q  column  2q+l 


D.  = 


Row  2f  -1 


Row  21 


^^2<+l^‘‘’^2f+2 


^^21 +2  ^■‘■^21  +  31 


f" 


^3f+2‘^2f +1 


^2l+4"^2f+2 


n-4 


n-3 


n-2 


• 

• 

•  • 

•  • 

(a  ^l+a  , 

'  n-2'  n-1 

^n"^n-2 

-(a  ,)+a 

n-r  n 

-^n-1  • 

(a  ) 

-a 

n 

n 

^2q-2f+l‘'’^2q+2f+l 

^2q-2i  +  2''’“2q+2f+2 

■*2q-2f-r^2q+2f  -1 

'^2q-2f  "^2q+2i 

^2q-2f '*’^2q+2/+2 

^2q-2f  +  l‘^^2q+2f  +  3 

■*2q-2f-2'®2q+2f 

•  . 

'*2q-2<-r^2q+2f+l 

It  should  be  noted  that  the  elements  in  column  q  and  row  i  of 

Dj,  are  identical  to  the  elements  in  column  q+1  and  row  f+1  of 

A  or  of  A  ,+B 

n-1  n-1 
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5.  We  subtract  row  n-2  from  row  n-3.  then  after  simplifying 
the  new  row  n.3.  we  subtract  it  from  row  n-4.  and  proceed  in  the 
same  fashion  up  to  the  first  row.  Then  we  add  column  1  of  the  new 
determinant  to  column  2,  then  after  simplification,  we  add  the  new 

column  2  to  column  3  and  repeat  the  same  process.  We  finally  obtain 
for  the  following: 


column  1  column  2  column  2q  column  2q+l 


row  1 

a  +a, 

0  3 

^2q-l'^^2q+2 

^2q^^2q+3 

row  2 

^4 

,+a 

2q-2  2q+3 

^2q-l'‘’^2q+4 

row2f  -J 

^2f +1 

^21+2 

^2q-2|  +  l'*'^2q+2f 

^2q-2|+2 

'^^2q+2f  +  l 

row  21 

^21+2 

^21+3 

^2q-2l+ 

^2q-2|+l 

— . 

. 

^2q+2l+l 

''■^2q+2£  +  2 

It  IS  noticed  that  the  new  expression  of  is  identical  to  the 

expression  of  A^.2'''®n-2’  shown  earlier  on  p.  (29).  Therefore 

^1  “  '^n-2^®n-2‘  ^^"'ilarly  we  can  show  the  identity  of  D  =  A 
-B  2  n-2 

n-2 

Therefore  with  this  equivalence  we  have  verified  the  following 
identity: 


or  equivalently  the  third  property 
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is  established. 

Discussion  of  the  first  property: 

The  first  property  as  given  in  Eqn.  25.  can  be  also  written  as 

■^k  “  '^k-l'^k+l  ■  ®k-l\+l  ’  =  2,  3,  4,  5. . .  n-1 


We  can  readily  verify  the  above  property  for  the  limiting  case.i.  e.  , 

when  k  =  n-1,  by  combining  the  second  and  third  properties  discussed 
earlier  as  follows; 

The  third  property  gives; 


•  •  *  V  '  '  >  ^„-2 


The  second  property  gives; 


A 

n 


(a  +a,+. 
o  2 


.  .  a_  +.  .  .  )(A  -B  ) 
2p  ''  n-1  n-r 


and 

i£  we  multiply  the  third  property  by  and  use  the  second 

property  we  obtain; 


n-1 


I  =  A  A  , 

ri”l  n  n-2 


B  B  _ 
n  n-2 


The  above  is  exactly  the  first  property  for  the  limiting  case.i.  e. , 
when  k  =  n-1.  By  actual  expansion,  the  first  property  has  also  been 
verified  for  k  =  2,  3,  4,  5.  In  order  to  complete  the  proof  it  has  to  be 
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shown  that  it  is  valid  for  any  k  between  5  and  n-1.  This  proof 
could  be  achieved  in  one  of  the  following  two  procedures: 

a)  By  determinant  manipulation  as  has  been  done  for  the  other 
properties,  if  we  write  for  the  first  property  the  following  e(]uivalent 
form: 


b)  By  induction  method,  i.  e.  ,  to  show  if  it  is  valid  for  k-1 ,  it  is  also 
valid  for  k. 

Both  the  above  procedures  involve  difficult  and  complicated 

manipulations  which  were  not  attempted  in  this  report.  However,  we 

may  present  a  simple  heuristic  argument  to  indicate  that  the  first 

property  holds  for  all  k.  This  is  based  on  the  following  observation. 

If  we  assume  a  certain  n,  i.  e.  ,  n  =  5,  then  the  stability  constraints 

are  given  as  1  >  0,  Aj<  0,  A^>  0,  A^<  0,  >0,  A^  <  0.  Now  if 

we  assume  any  general  n  >5,  the  stability  constraints  are  given  by 

1>  0,  A^<  0,  A^>  0,  A^<  0,  A^>  0,  A5<  0,  A^  >  0,  A^<  0... 

The  A's  for  the  general  case  up  to  A.  are  the  same  as  the  A's  for 

i) 

n  =  5  except  for  replacement  of  the  specific  n  =  5,  the  general  n  in 
computating  the  determinants.  Furthermore,  any  relationship  that 
holds  between  the  A's,  i.  e.  ,  A^-B^  =  ^-l^k+r^k-l^k+l* 
for  n=  5,  also  holds  for  any  n.  Thus  one  may  deduce  that.  "If 
the  first  property  is  verified  for  any  specific  n.  it  also  holds  for  any  n.  " 
Based  on  the  above  deduction,  we  can  use  the  limiting  case  of  the 
first  property  to  extend  the  range  for  n=6,  7.  .  .  For  instance  if  n=  7, 
then  the  first  property  verified  for  k=n-l,  becomes  also  valid  for  k=6 
and  for  all  n.  Similarly  we  may  proceed  step  by  step  in  the  same 
fashion  to  cover  all  the  intermediate  cases  of  k. 


*By  showing  that  this  equation  which  holds  for  k=5  and  n=6, 
to  be  valid  for  k=5  and  any  "n",  then  a  rigorous  proof  has  been 
constructed  by  using  the  limiting  case  to  extend  the  range  of  k. 
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Admittedly,  the  above  argument  doesn't  constitute  a  rigorous 
proof  but  only  indicates  a  convincing  argument  that  the  first  property 
cannot  be  violated  for  any  k  between  5  and  (n-1).  One  can  also 
use  the  expansion  method  to  verify  the  results  for  higher  "k",  however, 
this  again  involves  a  complicated  procedure. 

In  summary,  the  material  of  the  appendix  yields  rigorous 
proofs  for  the  second  and  the  third  property  and  from  these  pro¬ 
perties  a  heuristic  argument  for  the  validity  of  the  first  property 
is  indicated. 
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